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In this contribution. I discuss some basic techniques that can be used to simulate radiative 
transfer through porous media. As specific examples, I consider scattering transfer through a 
clumped slab, and X-ray emission line formation in a clumped wind. 



1 Introduction 



A porous medium is one that is macroclumped: the 
typical scale of the clumps is larger than the mean 
free path £ of photons. This represents the opposite 
end of the 'clumping spectrum' to a microclumped 
medium, for which the clump scale is smaller than 
£. The distinction between these two limits is 
more than academic; although microclumping usu- 
ally only affects the volume emissivity of a medium 
(which depends on the local density squared) , poros- 
ity/macroclumping can result in the medium having 
an effective opacity KcS that is significantly smaller 
than the microscopic value k of the material com- 
prising it. 

Recent investigations of porosity (e.g., Oskinova, 
Feldmeier & Hamann 2004; Owocki & Cohen 2006) 
have been focused primarily toward evaluating its 
observable consequences. As such, these studies 
have relied largely on very simple models for the ra- 
diative transfer. In this contribution, I discuss some 
of the basic techniques and codes that I have devised 
to move beyond these simplified models, allowing 
consideration of the full radiative transfer problem 
for a clumped medium. 

The standard approach to radiative transfer is to 
solve the radiative transfer equation (RTE), 



(1) 



where / is the specific intensity, S the source func- 
tion, and r the optical depth coordinate. However, 
solving the RTE is not the only way to tackle the 
problem. Instead of the /-based 'field' approach em- 
bodied by the RTE, one can adopt a 'particle' ap- 
proach, following individual photons as they prop- 
agate through a medium. This is precisely what 
Monte-Carlo radiative transfer entails, and in the 
following section I apply this approach to model 
scattering transfer through a clumped slab. 



2 A clumped scattering slab 

To investigate radiative transfer in a clumped scat- 
tering medium, I have developed freyr, a Monte- 
Carlo code for simulating the propagation of pho- 
tons through a plane-parallel slab composed of dis- 
crete, spherical clumps. The clumps arc randomly 
positioned, but each has the same radius Tc and a 
uniform density pc- Photons are introduced at the 
bottom of the slab with a random upward direction. 
A nominal 'optical distance to travel' r is assigned 
to each photon based on the expression 



\ogx, 



(2) 



where a; is a uniform random deviate in the range 
< a; < 1. This optical distance is converted into 
a physical distance i by inverting the optical depth 
equation 



(3) 



here, the microscopic opacity k is assumed to be 
constant, while the local density p varies with posi- 
tion as the photon passes in and out of spheres. Af- 
ter moving a distance i in the appropriate direction, 
the photon is scattered: a new random direction is 
chosen (corresponding to isotropic scattering), and 
a new optical distance is picked using eqn. (2). This 
procedure is repeated until the photon escapes from 
the top or the bottom of the slab; periodic boundary 
conditions are applied at the sides of the slab. 

Conceptually, this procedure is quite straightfor- 
ward; but from an implementation perspective the 
tricky part comes in inverting the optical depth 
equation (3). The key to this task lies in the recog- 
nition that along any given ray through the slab, the 
density is a piecewise-constant function. Changes in 
p arise only at those discrete points where the ray 
enters or exits a sphere. To locate these points, the 
parametric equation for the ray traversed by a pho- 
ton is written as 



r = ro + Mr. 



(4) 
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Figure 1: The ratio between the effective opacity Koff 
and the microscopic opacity k, plotted (di- 
amonds) as a function of clump optical 
thickness Tc- The solid line shows the cor- 
responding prediction of the simple bridg- 
ing law (14). 



Here, rg is the starting position vector of the pho- 
ton, and the unit vector dr gives the photon's di- 
rection. To establish whether this ray intersects a 
sphere located at rg, freyr calculates the following 
quantities: 



V = ro - 
B = 2dr • V 

C = vv - 



4C 



(5) 
(6) 
(7) 
(8) 



If £» > 0, 
locations 



then the ray pierces the sphere at the 



4 = 



(9) 



these correspond to the points where the density un- 
dergoes a jump. Between these points, p is constant, 
and T varies linearly - so it is straightforward to re- 
construct the function t{1) from eqn. (3), and then 
to invert this function to find £ given r. 

For each leg of a photon's journey, the ray-sphere 
intersection test (5-9) must be performed against 
every sphere contained in the slab. Potentially, this 
can be quite time consuming, and so freyr uses a 
spatial partitioning approach to reduce the number 
of tests that must be done. After the distribution 
of spheres is calculated (i.e., a set {y^} of random 
centers), the slab is divided up into a 3-dimensional 
array of voxels. Each voxel stores a index list of those 
spheres that fall part or all of the way inside it. To 
find out which spheres a ray intersects, it is then 



only necessary to test the subset of spheres that are 
indexed by the voxels that the ray passes through. 
Efficient algorithms exist (e.g., Amanatides & Woo 
1987) for quickly finding these voxels. 

Having described the basic principles of freyr, 
let me now present some illustrative results. Fig. 1 
plots the ratio Kbs/k between the effective and mi- 
croscopic opacities of a clumped slab, as a function 
of the averaged clump optical thickness 
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(10) 



To obtain the Koff datum at each abscissa, a 
FREYR simulation is used to determine an empiri- 
cal value for the slab transmittance 



(11) 



here, Ni^j is the number of photons injected at the 
bottom of the slab during the simulation, and A'osc is 
the number that eventually escape through the top. 
Formally, the transmittance is a function only of the 
slab optical thickness, so that 



(12) 



for some function /(). Thus, armed with the 
empirically-measured value of T, and the known val- 
ues of the slab mean density p and vertical extent 
Az, the effective opacity follows as 



pAz 



(13) 



The only difficulty is to determine the inverse of the 
transmittance function, f~^{), but this can be done 
numerically using Monte-Carlo simulations of a ho- 
mogeneous, isotropically scattering slab. 

Looking at Fig. 1, the onset of porosity is plain 
to see. Once the clumps become optically thick (i.e., 
Tc > 1), the effective opacity of the slab begins to fall 
well below the microscopic value n. This drop-off is 
well approximated by 



Kotf_ 
K 



C 



c 



(14) 



(shown in the figure by the solid line), where the con- 
stant C w 1.4. This expression is a modified form 
of the bridging law introduced by Cohen & Owocki 
(2006). In the limit of small Tc, it correctly repro- 
duces the microclumping case, for which KbB equals 
the microscopic opacity k. In the opposite limit of 
large Tc (i.e., macroclumping), the effective opacity 
has a scaling 



Tc 



3C 

4/0c?'c 



(15) 



Multiplying through by the clump mass Anpcrc/S 
then gives the effective cross section of the clumps. 
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Figure 2: X-ray line profiles for a r, = 3 clumped 
wind. The solid (dotted) profile corre- 
sponds to optically thin (optically thick) 
clumps; both profiles are normalized to 
have a maximum flux of unity. 



(7eff=C^r2. (16) 

This can be recognized as the geometric cross sec- 
tion 7rr^, scaled by the constant C. The basis for 
this result is that the clumps - when they become 
extremely optically thick - behave like individual 
scattering centers, with an interaction cross section 
equal to their geometric cross section. The appear- 
ance of the constant C is to correct for the angular 
scattering profile of each individual clump. If this 
profile were isotropic, then wc would obtain C = 1; 
but in fact the profile is that of a sphere with a 
Lambert-law surface (see Schoenberg 1929), leading 
to the 'observed' value C = 13/9 w 1.4. 

3 A clumped absorbing wind 

The utility of Monte-Carlo is that it works in cases 
such as scattering where the source function S is dif- 
ficult to obtain. If in fact S is already available, then 
a simple formal solution of the RTE (cf. eqn. 1) is 

j Sc-^dT. (17) 

Calculating the intensity / using this expression 
is invariably faster than performing an equivalent 
Monte-Carlo simulation. (This is a specific instance 
of a general rule: Monte-Carlo is guaranteed to be 
the slowest way to solve a radiative transfer prob- 
lem. Of course, sometimes - as with scattering - it 
is the only straightforward way.) 



The calculation of X-ray emission line profiles in 
clumped winds (see, e.g., the contribution by David 
Cohen) is a specific example of a case where the 
formal solution (17) can be used. X-ray line pho- 
tons arc emitted according to some simple paramet- 
ric function of radius, meaning that the source func- 
tion is easy to calculate. These photons are then 
re-absorbed by an opacity associated with a distri- 
bution of clumps. As with the Monte-Carlo simula- 
tions, the tricky part in modeling this re-absorption 
involves evaluating the optical depth r. 

Fortunately, for spherical clumps the same set of 
techniques (cf. §2) can be applied. If the inter- 
clump medium is a vacuum, then once again the 
density along any given ray is a pieccwisc-continuous 
function, and the ray-sphere intersection testing (cf. 
eqns. 5-9) can be used to reconstruct t(^). If, more 
realistically, there is a smooth inter-clump medium 
following a wind expansion law p = M / (Airr'^v), 
then the contribution of this medium to the optical 
depth along a ray can be added in using numerical 
quadrature. 

Fig. 2 shows some example X-ray line profiles cal- 
culated using BOREAS, a simple formal solver for a 
spherical-clump wind. For both profiles, the optical 
depth to the star (cf. Owocki & Cohen 2006) is 
3; however, the profiles differ in the nature of the 
clumping. The solid-line profile is for a configura- 
tion composed of many optically thin clumps, and it 
shows the characteristic asymmetric (blue-skewed) 
shape usually associated with a smooth wind. Con- 
versely, the dotted-line profile is for a configuration 
composed of fewer, optically thick clumps. With 
porosity acting to reduce the effective opacity of the 
wind, a much more symmetric line profile is seen in 
this latter case. 

4 Summary 

The spherical clumps considered herein admittedly 
rather idealized. However, it is precisely for such 
simple, idealized cases that one can hope to fully un- 
derstand the precise manner in which porosity mod- 
ifies the effective opacity of a medium. Both freyr 
and BOREAS - and the techniques they implement 
- will certainly prove useful in developing such an 
understanding. 
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